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I discuss the physical picture underlying the evolution equations with Pomeron loops recently 
derived in multicolor QCD at high energy and qualitatively explain the notion of 'self-duality'. 

The rapid rise of the gluon density in a hadronic wavefunction with increasing energy together 
with the phenomenon of gluon saturation open the way towards a perturbative treatment of 
high-energy scattering in QCD in the vicinity of the unitarity hmit ^ . To that aim, one needs 
evolution equations which include all the microscopic processes contributing in perturbative 
QCD to the growth of the gluon distribution and its eventual saturation. These processes have 
been discussed in a different article in these Proceedings'^ and are also summarized in Fig. 1. 
They include the BFKL evolution (Fig. l.b), which involves a 2 — > 2 gluon vertex and leads 
to the rapid growth of the gluon density with increasing energy, together with gluon splitting 
vertices 2 — > n with n > 2 (Fig. l.c) and gluon recombination vertices n — > 2 (Fig. l.d). One 
must emphasize here that all such processes — in fact, all the gluon number changing vertices 
m ^ n with arbitrary m, n > 2 — are important for computing scattering amplitudes at very 
high energy^, and this even when the projectile is a relatively simple object, like a color dipole 
(the interesting case for deep inelastic scattering at small x). Indeed, the non-linear effects 
responsible for gluon saturation and unitarization involve higher n~point correlation functions 
(cf. Fig. l.d), which in turn are generated through splitting processes in the low density regime 
(cf. Fig. l.c). In Fig. 1, we have also indicated some simple scattering processes involving one 
or two external dipoles that we shall use below to probe the gluon evolution in the target. 

Gluon-number changing vertices started to be computed in perturbative QCD in the early 
nineties^, but their systematic inclusion into evolution equations turned out to be a complicated 
problem, whose general solution is not yet known. The merging vertices are presently included 
in the functional JIMWLK equation ^, whereas the splitting vertices are rather described by 
the dual version of the JIMWLK Hamiltonian ^ (see below for the notion of 'duality'). But 
the complete evolution Hamiltonian, which should include both splitting and mergings, and 
therefore be self-dual^'^ , is not yet fully known (see however Refs. ^). 

The desiderata of including both splitting and merging has been nevertheless achieved ^' ^ 
within the multicolor limit Nc ^ oo, which allows for important technical simplifications and 
also for the use of a suggestive dipole language ^' (see also Refs. ^.n^is^^ 

Namely, at large Nc and in the dilute regime, the gluons in the target wavefunction can be 
effectively replaced by color dipoles, which interact with the external dipoles from the projectile 
via two gluon exchanges (compare in this respect the processes in Figs, l.a and 2. a). Then the 
BFKL evolution together with the 2 — > 4 gluon splitting (as previously depicted in Figs, l.b 

"At large Nc, and for a dipole projectile, one can ignore all the splitting vertices except for the vertex 2 — > 4. 




and c) can be represented as one target dipole splitting into two, followed by the interaction 
between the child dipoles and the projectile (cf. Figs. 2.b and c). Note that, whereas the 
scattering with a single external dipole serves as a measure of the average dipole number density 
in the target, and thus of the standard BFKL evolution (cf. Fig. 2.b), on the other hand, the 
scattering with two external dipoles (cf. Fig. 2.c) can measure both dipoles produced after a 
splitting, and thus probes the dipole number fluctuations in the target. 

But the dipole picture breaks down at high density^, and thus cannot be used to also 
describe the 4 — > 2 recombination process in Fig. l.d. However, by viewing the diagram in 
Fig. l.d upside down, it becomes clear that a merging process in the target can be alternatively 
interpreted as a splitting process in the projectile ; since the latter is dilute, the dipole picture can 
then be used to describe its evolution. The two diagrams in Figs. 3.b and c represent the same 
processes as previously shown in Figs, l.b and d, respectively, but now from the perspective of 
projectile evolution. In particular, the diagram in Fig. 3.c describes the simultaneous scattering 
of both child dipoles in the projectile, which is the unitarization mechanism at high energy. 
We see that the unitarity corrections correspond to either gluon saturation, cf. Fig. l.d, or to 
multiple scattering, cf. Fig. 3.c, according to the frame in which we choose to view the evolution. 

The previous considerations show that, in order to construct evolution equations which 
include both splitting and merging at large N^, it is in fact sufficient to know the elementary 
dipole splitting vertex, that is, the probability density for an original color dipole to split into 
two new dipoles per unit rapidity. This has been first computed by Al Mueller^, and reads 
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where dr is the rapidity increment (r ~ Ins is the rapidity gap between the projectile and the 
target), the pair {x, y) denotes the transverse coordinates of the quark and the antiquark in the 
original dipole, and {x,z) and {z,y) refer similarly to the two final dipoles. The latter have a 
common leg at z, which physically represents the transverse coordinate of the emitted gluon. 

By using Eq. (1) together with the picture of projectile evolution, cf. Fig. 3, one can imme- 
diately deduce the following evolution equation for the average scattering amplitude {T{x, y))r 
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Figure 3: Projectile evolution in the dipole picture. 

for a projectile made with a single dipole (a;, y) (with Os = asNc/ir) : 
d 



dr 



{T{x, y))r = M{x, y,z)(^- T{x, y) + T{x, z) + T{z, y) - {T{x, z)T{z, y)) ^. (2) 



Together, the first three terms in the r.h.s., which are linear in (T), represent the BFKL evolution, 
that is, the scattering of only one of the two child dipoles (cf. Fig. 3.b) and the reduction in the 
probability that the incoming dipole remains in its original state (the 'virtual term' —{T{x, y))). 
The last, non-linear, term {T(x, z)T(z, y))T describes the simultaneous scattering of both child 
dipoles, cf. Fig. 3.c. So long as the scattering is weak, (T)^ <C 1, this last term is negligible and 
the amplitude rises very fast (exponentially in r), according to the BFKL equation. But when 
(T)t- ~ 1, the non-linear term becomes important and tames the growth. Thus, as anticipated, 
the multiple scattering is the mechanism leading to unitarization. 

But in order to study this mechanism in detail, one also needs the equation satisfied by the 
2~dipole amplitude {T{x, z)T{z,y))T-, or, more generally, {T{xi,yi)T{x2,y2))T- This involves, 
first, an evolution similar to Eq. (2) for any of the two incoming dipoles {xi,yi) or (0^2,1/2) 
(while the other one is simply a spectator). This will contribute terms linear in {TT)r (the 
BFKL terms) together with unitarity corrections proportional to (TTT),-. Besides, the two 
external dipoles can 'see' a fluctuation in the gluon in the target, that is, they can directly probe 
the 2 — > 4 splitting vertex in Fig. I.e. To estimate this additional effect, it is preferable to adopt 
the point of view of target evolution, cf. Fig. 2.c, which implies^ 
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The r.h.s. of this equation should be read as follows: A dipole {u, v) from the target splits into 
two new dipoles {u,z) and {z,v) with probability {as/2iT)A4{u,v, z) (cf. Eq. (1)), then the 
two child dipoles scatter off the external dipoles, with an amplitude cks^o for each scattering. 
Finally, the 'amputated' amplitude (l/ag)V^V^ {T{u,v))^ is, up to a normalization factor, the 
average dipole number density {n{u,v))^ in the target. 

Note that, although suppressed by an explicit factor a^, the contribution of the dipole 
number fluctuations to the evolution of {TT)r is in fact a dominant effect in the very dilute regime 
where (T),- ^ : Indeed, in this regime {T){T) < a^{T), hence the dominant contribution to 
the rise of {TT)r comes from fluctuations, via Eq. (3). This conclusion is in agreement with the 
qualitative analysis of the fluctuations in the companion article, Ref. . 

It is now clear that the evolution equations for dipole scattering amplitudes form an hierar- 
chy, in which the general equation reads as follows (with t = a^r and T*^""-* = (T(1)T(2) . . . T{n))T-) 



^a'iAo{x,y\u,v) is the standard amplitude for dipole-dipole scattering via two-gluon exchange" 



in schematic notations which ignore combinatoric factors and the non— locahty of the various 
vertices in the transverse space. The first two terms in the r.h.s. correspond to the BFKL 
evolution and the unitarity corrections, respectively, while the last term describes dipole number 
fluctuations in the target wavefunction. The iterative solution to these equations can be given a 
diagrammatic representation in terms of 'BFKL Pomerons' (the Green's function of the BFKL 
equation) which split and merge with each other, and thus form Pomeron loops^. 

An hierarchy like that in Eq. (4) can be exactly represented by a Langevin equation with 
a specific noise term. In statistical physics, this equation is known as the stochastic FKPP 
equation (sFKPP), and has the following generic structure (see Ref. ^ for details) : 



^ = T-T'^ + ^fi^v{t) with {v{t)v{t'))=5{t-t'). (5) 

Such an equation represents a convenient starting point for numerical simulations, and particular 
limits of it have been already studied in this way . 

We shall conclude this brief presentation with a qualitative explanation of the concept of self- 
duality in the context of the high-energy QCD evolution 5-'^. The seemingly trivial fact that a 
same n — > 2 vertex can be viewed as either merging in the target or splitting in the projectile has 
the physically important consequence that, in order to be hoost invariant^ the evolution equations 
for the scattering amplitudes must include the effects of both splitting and merging. Indeed, 
physical quantities must come out the same whatever is the system that we choose the evolve 
(the projectile or the target), since such a choice is tantamount to using a particular Lorentz 
frame. Let's then assume that the evolution Hamiltonian H — the operator which governs 
the one-step evolution of the gluon fields in a hadronic wavefunction — contains a particular 
merging vertex n ^ 2. When acting on the gluon fields in the target, this particular term in 
H generates an evolution for the scattering amplitudes which is equivalent to that that would 
be produced by a specific splitting vertex (the mirror, or 'dual', image of the original merging 
vertex) acting on the gluons in the projectile. Hence, in order to yield the same physical result 
when directly applied to the projectile, the Hamiltonian must also include the specific splitting 
vertex mentioned above, and therefore be self-dual. The precise duality transformation relating 
merging to splitting vertices depends upon the detailed structure of the scattering operator, and 
can be found in Refs.^~^. The Hamiltonian underlying the evolution equations with Pomeron 
loops presented above is indeed self-dual^, and so is also the more general Hamiltonian in Refs.^. 
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